SPECIFICATION 



Electronic Version 1.2.8 
Stylesheet Version 1 .0 

SCATTER CORRECTION METHOD 
FOR NON-STATIONARY X-RAY 
ACQUISITIONS 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH & 
DEVELOPMENT 



g [0001] The government may have rights in this invention pursuant to Subcontract MDA 

905 001 0041 issued from the Office of Naval Research/Henry M.Jackson Foundation. 



m 



BACKGROUND OF THE INVENTION 



[0002] The present invention relates generally to scatter-correction for X-ray or the like 

fjj type of examination technique, and more specifically to a scatter-correction technique 

iff I 

M ? for situations wherein the angle of incidence of an incident ray from a radiation source 

III 

m impinges on a detector at an angle which is not exactly zero degrees or exactly at 

fM right angles to the detector, and thus allows for the source to be moved with respect 
to the detector or vice versa. 

[0003] Breast cancer is a significant cause of mortality among women in the United 
States. Currently the most widespread breast-cancer screening technique is X-ray 
mammography, where malignancies are detected based on their characteristic shape 
(typically speculated) and radiological density (similar to breast glandular tissue). 

[0004] Early detection of breast cancer generally improves a patient's chance of survival 
and is thus of critical importance. The task of detection of malignancies is 
complicated by the fact that the breast is a radiologically complex 3D structure that is 
imaged by 2D X-ray projection. Thus, malignant lesions can be masked by 
surrounding glandular tissue, decreasing contrast and distorting the lesion's perceived 
shape. A significant number of malignancies are missed in mammography because of 
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the presence of this increase of structural noise brought about by the projection of 3D 
structures in the breast onto a 2D image. 



[0005] Tomosynthesis is a promising alternative to conventional two-dimensional 

mammography. This imaging technique resides in the acquisition of plural low-dose 
2D X-ray projection images, where the X-ray source is moved along an arcuate 
trajectory with respect to the stationary breast and/or detector. X-ray projections are 
thus acquired at different angles. These angular projections are then combined in a 
3D image reconstruction of the breast. In the reconstructed image, the overlap of 3D 
structures is mitigated, contrast is enhanced and the morphological features of 
malignancies (spiculation) can be better preserved than in 2D X-ray imaging, 
therefore leading to improved sensitivity and specificity with respect to 2D X-ray 
*■* mammography. 



CO 

o 
m 



f| [0006] In tomosynthesis, the use of a non-stationary X-ray source precludes the use of a 
traditional grid for scatter rejection. With larger source angles with respect to the 
detector normal, the scatter signal will increase due to the longer distances that the 
X-rays traverse through the breast tissue. Therefore, there is a need to be able to 

h4 estimate the scatter component on an X-ray image that has been acquired and 

s t% 

fa compensate for the scattering in order to improve the resolution/effectiveness of the 

in 

|L| examination of the breast tissue. 

W " BRIEF SUMMARY OF THE INVENTION 

[0007] More specifically, a first aspect of the present invention resides in a method of 

estimating scatter in an image, where radiation from a radiation source is transmitted 
through an object to be imaged and is incident on a detector at an inclined or "non- 
zero" angle. This method comprises explicitly modeling an effect of angular incidence 
of the radiation on a scatter signal, including deriving, using an empirically based 
correction which includes the effect of a non-zero angle of incidence of the radiation 
on the detector, an estimate of the scatter signal in image data derived from the 
incident radiation on the detector. 



[0008] 



A second aspect of the invention resides in a method of calculating and removing 
scatter from an image, where radiation emitted from a source is transmitted through 
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an imaged object, is incident on a detector at a non-zero (acute) angle, comprising: 
explicitly modeling the effect of angular incidence of radiation in a scatter signal; re- 
normalizing the scatter signal depending on the angle of incidence of radiation and 
the thickness of the object imaged; and correction of scatter from an image based on 
inverse-filtering. 

[0009] A third aspect of the invention resides in a method of estimating scatter in an X- 
ray image, where radiation emitted from an X-ray source is transmitted through an 
imaged object and is incident on the X-ray detector at an acute or non-zero angle, 
comprising: explicitly modeling the effect of angular incidence of X-ray radiation in a 
scatter signal; and calculating, via convolution, an estimate of the scatter signal in an 
X-ray image. 

Q [0010] A fourth aspect of the invention resides in a method of estimating an 



tit asymmetrical scatter kernel wherein asymmetry is introduced by angular incidence of 

m radiation, which has been emitted from a source and transmitted through an object to 

St 

f§ be imaged, on a detector, comprising: modifying scatter that would be derived 



m 



wherein the radiation is directly incident on the detector with zero degrees of 



Q inclination, using an asymmetry factor which indicates the shape and magnitude of 

ifi ■ 
the scatter kernel and which varies with an angle at which the radiation is incident on 



the detector, a mean attenuation coefficient of the object, and a distance the radiation 



m ' 

?.:r v 
B 

fU has traveled through the object. 

Hi 

[001 1] Further aspects of the invention respectively reside in computer readable media 
which are encoded with programs executable by a computer for implementing the 
above mentioned methods. Additionally imaging systems which include hardware and 
computer software and which are adapted to implement the above mentioned 
methods, also comprise further aspect of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[001 2] Figure 1 is a schematic representation of an example of an X-ray system that the 
embodiments of the present invention can be applied. 

[001 3] Figure 2 is a diagram which depicts the foundation of equations 1-14. 
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[0014] 



Figure 3 is a schematic perspective view showing an experimental setup for 
angular measurements of scattered radiation. 



[001 5] Figures 4-9 are graphical comparison between open-geometry (y ), scatter- 

K 

free (b ) and scatter-corrected (b f ) signals for a 25 ° angular acquisition. 

K. K 

[001 6] Figures 1 0 and 1 1 are histograms wherein the histogram of Fig. 1 0 shows the 

percent deviation between the scatter-free signal and the open geometry signal within 
the 5 cm breast phantom (25 0 acquisition), and Fig. 1 1 shows, after scatter 
correction, the percent deviation between the scatter-free signal and the scatter 
corrected signal within the 5 cm breast phantom (25 ° acquisition). 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

p [001 7] The present invention is directed to a scatter-correction technique which enables 
5 the angle of the incident radiation to be taken into account and thus includes the 

Ul relative position between the radiation source and the radiation detector arrangement 



m 



which receives the radiation in the calculation and correction for scatter. It will be 



noted that the angle of incidence can very closely approach zero degrees (viz., exactly 



g parallel to the normal of the detector plane and thus effectively encompass all angles 

iff 

which are encountered and in which scattering occurs and needs to be either 
1ft determined and/or corrected for. 

n$ [001 8] Before dealing with the concept on which the present invention is based, it is 

appropriate to consider an example of the type of device to which the concept of the 
invention can be applied. 

[0019] Figure 1 is a schematic representation of an imaging system 10, such as 

conventionally used in X-ray mammography imaging and tomosynthesis imaging. This 
system 10 includes at least one radiation source 14, such as an X-ray source, and at 
least one detector array 16 for receiving radiation emitted from the radiation source at 
one or more projection angles. This system can be adapted, for example, so that the 
radiation source 14 projects a cone-shaped beam of X-rays which pass through 
object 1 2 and impinge on the detector array 1 6. 



[0020] 



The detector array 1 6 has a plurality of detector elements (not shown per se) 
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arranged in rows and columns such that image data of the object 12 of interest can be 
suitably collected. As is conventional the detector elements can each comprise one or 
more photosensors, such as photodiodes. The radiation which is incident on a 
scintillator material enables the photosensors to measure, byway of change in the 
charge across the diode(s), an amount of light generated by X-ray interaction with the 
scintillator. 

[0021] Simply byway of explanation, the operation of radiation source 14 is governed by 
a control mechanism 28 of the imaging system 10. This control mechanism 28 
includes a radiation controller 30 that provides power and timing signals to the 
radiation source 14 and a motor controller 32 that controls a respective movement 
speed and position of radiation source 14 and detector array 16. 

y [0022] A data acquisition system (DAS) 34 which is included in the control mechanism 28 

P samples digital data from detector 1 6 for subsequent processing. An image processor 

hi 

VI 36 receives a sampled and digitized projection data set from DAS 34 and performs 

Ul 

Q image processing and scatter correction, as described herein. A system controller 38, 

m 

such as a computer or a microprocessor control the other elements of the system 1 0. 

Is??.:? 

if! [0023] Figure 2 is a diagram which demonstrates the various parameters in which 

i • 

|s.| following equations are couched. Figure 3, on the other hand, is a schematic 

\4 perspective view of an experimental arrangement which was used to conduct tests to 



ry 



verify the voracity of the effect of the algorithm that can implement the estimations 
enabled using the techniques outlined below. In this arrangement, a plurality of acrylic 
plates are disposed one on top of the other. 

[0024] This arrangement radiologicially approximates tissue such as breast tissue and 
thus enables a mammographic simulation. Although not specifically indicated, the 
arrangement is intended to represent, merely by way of example, a compressed breast 
thickness, which is included along with a detector cover thickness and a compression 
paddle thickness. This enables a thickness value to be determined and used in 
calculations to determine a length of the path along which a ray has passed and which 
length is used in the disclosed embodiment of the scatter correction process. 

[0025] Scatter correction model 
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[0026] A mammography image can be thought of as formed by unattenuated X-ray 
photons arriving directly from the source (direct events), plus photons that scatter 
into the detector from the breast and the acquisition system. The measured image y is 
then defined as 

y=b+s 9 Equation 1 

where b represents the direct events detected and s is the image formed by x-ray 
photons scattered. Here, we are interested in recovering the direct image b. The term 
s can be modeled as a convolution between a scatter kernel p and the direct image b, 
so that 

s = p*b y Equation 2 

where * is the convolution operator. The scatter kernel p can be though of as the 
y, probability that an X-ray photon directed onto pixel k of the detector will actually be 

:!:f scattered into pixel j. 

& 

y [0027] From Equations 1 and 2 we have 

m 



y=b+p*b. Equation 3 

n [0028] Explicitly, for any given pixel j in the experimental image, we can express 

a 

Equation 3 as 

N 

yj= b i + H Pith* 

where the index k spans all pixels of the detector (1 through N). 



m 

in 

.«? i; 



t i '- 
■ ■■■■■ 



[0029] 



Assuming that p is known for a given image acquisition system configuration and 
object being imaged, we can make an initial estimation of the scatter signal through 
the expression 

(o) P*y p*(b + p*b) p*b*(S + p) _ . . 

s { } = f— r = i — ~ — - — ' = — — , Equation 4 

1+^ \ + p i + p 

where 5 denotes the Kronecker delta function and ^ is the scatter kernel norm. 

Rearranging Equation 1 we can write an initial estimate of the direct events 

b^^y-s™. Equations 

The scatter corrected image can then be estimated iteratively using the expression 

=£+(Z^s-£)*i,^y, Equation 6 
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where we must set the scatter kernel norm to P<^ in order to insure convergence. 
[0030] Scatter Kernel Asymmetry 

[0031] The angular incidence of the X-ray photons introduces an asymmetry in the 

original scatter kernel Pjk . Originally, if the angle of incidence between the X-ray and 
the normal to the detector surface is zero degrees, the probability of a photon from 

ray k to scatter and be detected in pixel j, is assumed symmetric and spatially 

invariant with respect to the indices k and j, so the kernel could be written explicitly 
as Pj " k . Therefore the kernel p depends only on the relative distances between pixels k 
and j. This symmetry is readily seen in Figure 2, where the kernel Pjk ~ p s* . However, if 
the angle of incidence of the radiation is a *° , then the length of the photon path 
from the point of scatter is different for pixels j and j\ Consequently the probability of 
transmission of the X-ray photons scattered into pixels j and j' will differ and in 
general p^p h . 



I [0032] We introduce a correction factor which relates the scatter kernels at zero degrees 
W with the angular scatter kernel so that 



p JJt & N0 jk p Jjt , Equation 7 

wherein: 



Pjk is the shape and magnitude of the scatter at zero degrees; 



m 

III N is an empirically determined normalization factor; and 



p I . &> A is an asymmetry factor given by 

0 - 1 t Equation 8 

rz 

[0033] where m is the mean attenuation coefficient of X-ray photons in the breast tissue. 
[0034] Distance r is given by 



r = ., Equation 9 

where p the distance between pixels k and j. 

[0035] The distance r' is given by 



[0036] 



r = yl(p-z twaf +z 2 , Equation 1 0 

The estimated height z* (angular incidence) at which the scatter events take place 
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are given by 



where 



z = (/* - - M )cos(a) , Equation 1 1 



£ = — Equation 1 2 

cos(a) 



5SBIV 

E .. "8 



and y is the attenuation coefficient for the mean energy of the X-ray spectrum in 
acrylic (viz., the material used in the plates shown in Fig. 2 to emulate breast tissue). 

[0037] The estimated height z (perpendicular distance) can be obtained from the 
expression for z f by setting <x =0. 

[0038] The normalization factor N is given by 

N = Equation 13 



where f(t, oc ) is an empirically-determined scatter-kernel normalization factor that 



y depends on the thickness of the object imaged and the angle of the radiation source 

iff 

with respect to the normal of the detector plane. Equation 6 then remains 



III b ^ H ^ S -p* b < 



Equation 14 



III Determination of the scatter-kernel normalization factor 

|3 

[0039] The angular incidence of radiation introduces an asymmetry in the scatter kernel, 
as explained in the previous section. In addition, the magnitude of the kernel changed 
in a manner that was unaccounted for by the simple increase of path length 4 of 
photons through the imaged object as a function of <x . This introduces an empirical 
normalization function f(t, a ), which is determined by a series of measurements 
spanning object thicknesses t from 1 to 8 cm and radiation incidence angles from 0 0 
to 45 0 . 

[0040] Two images were taken for any given object thickness t and incident X-ray angle 
oc : The slit collimator image was taken as shown in Figure 3. The slit collimator was 
made of a sheet of Tungsten with a thickness of 20 mils. A 1 mm wide slit allowed X~ 
rays through. The signal directly beneath the slit is effectively scatter-free. 
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[0041] 



An open -geometry image was acquired in the same way as the slit collimator 
image, except that the collimator was removed, thus allowing for scatter in the image. 



[0042] The scatter-free signal b acquired in the slit-collimator image was then 

K 

compared to the scatter-corrected open-geometry signal b 1 by calculating the 

K 

reduced chi-squared x 1 between the scatter-free and scatter-corrected signals; 
l N 

X 2 ■ "^> 2 ' Equation 15 

An empirical kernel normalization factor was determined by minimizing for x 2 a given 
t and (X . A surface function f(t, a ) was then generated by regression. 

[0043] Validation through anthropomorphic phantoms 

[0044] Experimental data were acquired using a prototype tomosynthesis X-ray system 
y ;; for mammography. Three anthropomorphic breast phantoms were used for validation 

of the scatter-correction technique, with nominal thickness of 4, 5 and 6 cm. Their 
ff! glandular/fat percent composition was 50/50, 70/30 and 80/20, respectively. A full 

Q 

m 



set of projections (1 1 images) were acquired for slit collimator and open geometry for 
each phantom. Each set of projections covered tube angles ranging from +25 ° to -25 
° in 5 ° steps. The technique used was Rh/Rh, 32 kVp and 1 0 mAs per projection 
image. 

m 

'in- [0045] The open-geometry projections were scatter-corrected using the new scatter 
£1 kernel P jk , and compared to the scatter-free signal acquired using the slit collimator. 

m ' 

[0046] Figures. 3-8 show the open-geometry, scatter-free and scatter corrected signals 

for pixels directly below collimator slit. As will be appreciated, these figures show a 

comparison between open-geometry (y , ), scatter-free (b , ) and scatter-corrected 

k k 

(b* ) signals for a 25 ° angular acquisition. The mean difference between the scatter- 
free and scatter-corrected is less than ~ 3%. 



[0047] 



In the histogram shown in Figure 1 0, the % deviation between the scatter-free 
signal and the open geometry signal within the 5 cm breast phantom (25 ° acquisition) 
is shown. After scatter correction the percent deviation between the scatter-free 
signal and the scatter corrected signal within the 5 cm breast phantom (25 0 
acquisition) is show in the histogram of Fig. 10. 
[tl] 
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Table 1 



25 deg 


20 


deg 


15 dee 




scatter 


no scatter 


scatter 


no scatter 


scatter 


no scatter 




correction 


correction 


correction 


correction 


correction 


correction 


% bias 


o.ai 


-68.64 


3.38 


-67.16 


2.40 


-67.26 


^standard 


4.67 


16.90 


4.74 


15.40 


4.59 


14.06 


deviation 













10 


deg 


5 deg 


0 


deg 


scatter 


no scatter 


scatter 


no scatter 


scatter 


no scatter 


correction 


correction 


correction 


correction 


correction 


correction 


1.97 


-66.63 


1^9 


-65.86 


1.90 


-65,40 


4.47 


13.40 


4,57 


t.3,70 


4,59 


14.26 



m 



m 

If! 
m 

m 

m 



[0048] As will be appreciated, the scatter-corrected projections obtained with the present 
invention are comparable with the nearly scatter-free signal obtained using a slit 
collimator. Experimental data shows that the mean percent difference between the 
scatter corrected and the scatter-free projection signal is less than 3% in all cases, 
even in the edge region where there is rapid change in the thickness of the breast 
phantoms. 

[0049] As will be appreciated, this table shows the bias and standard deviation of the 

scatter-free signal and the open-geometry signal with and without scatter correction 
and highlights the effectiveness of the technique. 

[0050] The content of an article entitled Digital Tomosynthesis in Breast Imaging by 

LorenT. Niklason et al., Radiology. 205(2), 1997 p 399-406, is hereby incorporated in 
its entirety by reference thereto. 

[0051] The content of copending United States Patent Application No. 10/062,338 
(Attorney Docket No. RD 28622) entitled METHODS AND APPARATUS FOR 
CORRECTING SCATTER, is also incorporated by reference thereto. 



[0052] 



While the invention has been described in connection with a tomographical 
examination using X-rays, it will be understood that the estimation of scatter concept 
on which the present invention is based, is not limited to the use of X-rays nor living 
tissue and can be applied to inanimate objects such as air-craft wings or the like 
carried out using any suitable radiation source such as a gamma or electron beam 
source. Further, the various modifications and changes to the embodiment of the 
invention disclosed will be self evident to the those skilled in the art to which the 
present invention pertains or mostly closely pertains, when equipped with the 
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preceding disclosure and the need to modify and/or adapt the disclosed concept. 
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